arXiv:1507.04656vl [astro-ph.HE] 16Jul2015 


Compact Stars in the QCD Phase Diagram IV (CSQCD IV) 
September 26-30, 2014, Prerow, Germany 

http://www.ift.uni.wroc.pl/~csqcdiv 


Proton gaps and cooling of neutron stars with a 

stiff hadronic EoS 


Hovik Grigorian}'"^, Dmitry N. Voskresenski^ and David Blaschke^'^ 

^Laboratory for Information Technologies, Joint Institute for Nuclear Reseaarch, 
Joliot-Curie Street 6, If 1980 Dubna, Russia 

^Department of Physics, Yerevan State University, Alek Manukyan Street 1, 

Yerevan 0025, Armenia 

^National Research Nuclear University (MEPhI), Kashirskoe Shosse 31, llSfOO 
Moscow, Russia 

^Bogoliubov Laboratory for Theoretical Physics, Joint Institute for Nuclear Research, 
Joliot-Curie Street 6, If 1980 Dubna, Russia 

^Institute for Theoretical Physics, University of Wroclaw, Max Born Place 9, 50-20f 
Wroclaw, Poland 


Abstract 

The recent measurements of the masses of the pulsar J00737-3039B and of the 
companion J1756-2251 and pulsars PSR J1614-2230, PSR J0348-0432 demonstrate 
the existence of compact stars with masses in a broad range from 1.2 to 2 Mq. 
To fulhll the constraint Mmax > 2 Mq and to demonstrate the possibility of cooling 
scenarios for purely hadronic and further for hybrid stars we exploit the stiff DD2 
hadronic equation of state producing a maximum neutron star mass M ~ 2.43M0. 
We show that the ’’nuclear medium cooling” scenario for neutron stars comfortably 
explains the whole set of cooling curves just by a variation of the star masses without 
the necessity for the occurrence of the direct Urea reaction. To describe the cooling 
data with the very stiff DD2 equation of state we select a proton gap prohle from 
those exploited in the literature and allow for a variation of the effective pion gap 
controlling the efficiency of the medium modihed Urea process. Fast cooling of young 
neutron stars like it is seen in the data for Gas A is explained with the DD2 equation of 
state when the following conditions are provided: the presence of an efficient medium 
modihed Urea process, and a large proton gap at densities n L 2no vanishing for 
n L (2.5 — 3)no, where no is the saturation nuclear density. 
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1 Introduction 


Experimental data on surface temperatures of neutron stars (NSs) provide us with 
information about the neutrino emissivities of various processes, depending on the 
density behaviour of the NN interaction amplitude, values and density prohles of 
the proton and neutron paring gaps, the heat transport and the equation of state 
(EoS) of NS matter. Recently the situation has been improved with the observation 
of a segment of the cooling curve for the young NS in the remnant of the historical 
supernova Cassiopeia A (Cas A) [HE], with known age (~ 330 yr), for which the 
temperature and the rate of cooling have been followed over the past 13 years since 
its discovery [SlIllslEllIllHI. The data require the existence of a fast cooling process 
in the Cas A NS interior. On the other hand, the NS cooling model must also explain 
the compact object XMMU J173203.3-344518 [5] in a supernova remnant, for which 
the surface temperature has recently been measured. This object is hotter and older 
than Cas A, at an age between 10 and 40 kyr. Moreover, there exists information on 
surface temperatures of many other NS sources. It is not easy to appropriately explain 
these essentially different surface temperatures of various objects within the so called 
“minimal cooling paradigm”, where the only relevant rapid process is the so called 
pair-breaking-formation (PBF) process on 3 P 2 paired neutrons. The solution of the 
puzzle might be associated with a strong medium dependence of cooling inputs, as 
provided by the density (NS mass) dependent medium modihcations of the nucleon- 
nucleon interaction caused by the softening of the pion exchange contribution with 
increase of the density, and by the density dependent superfluid pairing gaps, see 
Hnum for details. The key idea formulated long ago [12] is that the cooling of 
various sources should be essentially different due to the difference in their masses. 
At that time there was the opinion that all NS masses should be hxed closely to 
the value of 1.4M0. The recent measurements of the masses of the pulsars PSR 
J1614-2230 [13], PSR J0348-0432 [TT] and J00737-3039B [13] and of the companion 
of J1756-2251 [in] have provided the proof for the existence of NS masses varying in 
a broad range (at least from 1.2 to 2 Mq). 

The most efficient processes within the “nuclear medium cooling” scenario are the 
medium modihed Urea (MMU) processes, e.g., nn —)■ npeu [121 [H], and the PBF 
processes N —)■ NpairVV, N = n or p [TH [19]. While being enhanced owing to their 
one-nucleon nature DSlEn], the latter processes are allowed only in the presence of 
nucleon pairing. Based on the assumption that the mass distribution of those objects 
for which surface temperatures are measured is similar to the one extracted, e.g., from 
a population synthesis, the very efficient direct Urea reaction, n —)■ npeu, should be 
forbidden in the majority of the former NSs, see [211E2]- The influence of in-medium 
effects on the NS cooling was hrst demonstrated in [23] with various EoS. The nuclear 
medium cooling scenario which was systematically developed further in [ 2111231123 ] 
provides a successful description of the whole set of known cooling data for NSs with 
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low magnetic fields. An overall fit of the data is obtained in our model for a strongly 
suppressed value of the 3 P 2 neutron pairing gap, thus being in favour of results by 

m- 

Reference [2nj has demonstrated an appropriate £t of the Cas A cooling curves 
with results from ACIS-S instrument yielding a surface temperature decline of 3 ... 4% 
over 10 yrs. For that description the lepton heat conductivity has been suppressed 
artihcially by a factor ~ 0.3 in the most favourable case compared to the result of 
Ref. [2H], which has been exploited in our previous works. The nucleon contribution 
to the heat conductivity is suppressed in our case by medium effects compared to that 
used in Ref. [2Hj. In the more recent work m we have used the result of Ref. [30] 
for the lepton heat conductivity which includes Landau damping effects. We have 
demonstrated that an appropriate fit of the Cas A data was then possible without ap¬ 
plying any artihcial suppression of the lepton heat conductivity. The HHJ EoS which 
has been exploited in our previous works has then been stiffened in [2n| for n > 4no to 
comply with the constraint that the EoS should allow for a maximum NS mass above 
the value M = 2.01 ± 0.04 Mq measured for PSR J0348-I-0432 by [Tl|, see also [T3] . 
However, the resulting EoS (labeled as HDD) which produces Mmax = 2.06 Mq might 
still not be sufhciently stiff, since the existence of even more massive objects than 
those known by now [131 El is not excluded. Incorporating systematic light-curve 
differences the authors of 1311 have estimated that the mass of the black-widow pulsar 
PSR J1311-3430 should at least be M > 2 .IM 0 . Furthermore, a deconfinement tran¬ 
sition in the NS interior would contradict these NS mass measurements, if one would 
use a soft hadronic EoS, cf. [32l[33l[33|. Therefore, the investigation of the possibility 
of hybrid stars requires a stiff hadronic EoS. Also, recent radius determinations from 
timing residuals suggest larger radii (albeit still with large uncertainties) and thus 
motivate the usage of a stiffer hadronic EoS [35] . 

Note that the authors of the recent paper [36] have explained the Cas A ACIS-S 
data for the NS mass M = 1.44Mq within the minimal cooling scenario using the 
BSk21 EoS, a large proton gap and a moderate 3 P 2 neutron gap. Hottest and coldest 
objects, however, can hardly be explained appropriately within the same scenario. 
Note also for completeness that the authors of Ref. 1371 have argued that the decline 
extracted from the ACIS-S graded mode data might be too steep. In spite of this 
uncertainty in the analysis of the data we, as well as Ref. [37], focus below mostly on 
a comparison with the ACIS-S graded mode data for Cas A. 

In the previous contribution [33] we have demonstrated preliminary calculations of 
the cooling curves within our nuclear medium cooling scenario exploiting the stiff DD2 
EoS [ 39 ] . However, in that work we have not performed any tuning of the parameters 
of the hadronic model. Additionally, as an alternative to the purely hadronic scenario, 
we have incorporated in [33] the possibility of a deconhnement phase transition from 
such a stiff nuclear matter EoS in the outer core to color superconducting quark 
matter in the inner core. Here we continue our study within the hadronic scenario 
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exploiting the stiff DD2 EoS by tuning the proton pairing gaps and the effective pion 
gap in order to construct a better fit of the cooling data. The deconfinement phase 
transition will be considered elsewhere. 


2 EoS of hadronic matter 

In our previous works [211 [211 US] we have exploited the HHJ (5 = 0.2) fit [ID] of the 
APR EoS jH]. While the latter EoS produces an appropriate maximum NS mass 
of M = 2.2 Mq, the HHJ EoS £t introduces an additional parametric correction of 
the high-density behaviour in order to avoid a causality breach. This comes at the 
price of a lowering of the maximum NS mass to M = 1.94 Mq, below the values for 
the measured masses of the pulsars PSR J1614-2230 [13] and PSR J0348-0432 [Tl] . 
Recently we have modihed this EoS by invoking an excluded volume for nucleons 
|29j . The so constructed HDD EoS stiffens for higher baryon density n resulting in 
an increase of the maximum NS mass up to the value Mmax = 2.06 Mq. 

However, a recent radius determination for the nearest millisecond pulsar PSR 
J0437-4715 [35] indicates that a still stiffer EoS might be needed to support (at 
2cr confidence) radii > 13 km in the mass segment between 1.5 and 1.8 Mq. The 
density dependent relativistic mean-field EoS of Ref. [12] with the well calibrated 
DD2 parametrization [39] meets this requirement. Moreover, the DD2 EoS fulfils 
standard constraints for symmetric nuclear matter around the saturation density 
and from nuclear structure. The density dependent symmetry energy agrees with 
the constraint by Danielewicz and Lee [H] and with ab-initio calculations for pure 
neutron matter [H]. The direct Urea reaction threshold is not reached within the 
DD2 EoS. However, due to the stiffness of DD2 EoS, it does not fulhl the ’’flow 
constraint” 051 for densities above 2no. This is the price to be paid for the possibility 
to increase the maximum mass and the radii of NS. 

The dependences of NS masses on central baryon density Ucen are shown in the 
left panel of Fig. [^for the HDD and DD2 hadronic EoS. We see that for a fixed NS 
mass the stiffening of the EoS leads to a redistribution of the density prohle in the NS 
interior so that the central densities get lowered. As a consequence, a slower cooling 
is expected for stars with the DD2 EoS when compared with stars of the same mass 
described by the HDD EoS. In other words, in order to cover the same set of cooling 
data with a stiffer EoS the range of masses attributed to the set of cooling curves shall 
be shifted to higher values. The right panel of Fig. [T] demonstrates the mass-radius 
relation (crust is not included). The stiffer DD2 EoS produces a larger radius than 
the softer HDD EoS. 

We shall now discuss the results for NS cooling obtained with both, the HDD and 
the DD2 EoS. 
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Figure 1: Mass vs. central baryon density (left) and vs. radius (right) for the stiff 
DD2 hadronic EoS (dash-dotted lines) and for the HDD EoS (solid lines). Crust is 
not included. 

3 Cooling 

In our previous works |21l EH ESI EH] we have demonstrated that the cooling history 
is sensitive to the efficiency of the medium modihed Urea process controlled by the 
density dependence of the effective pion gap shown in Fig. 1 of |2l] and to the value 
and the density dependence of the IS'o pp pairing gap. The results are very sensitive 
to the 3 P 2 nn gap. In our works we follow the analysis of ^7\ where this gap turns 
out to be negligibly small. Our results are rather insensitive to the treatment of the 
IS'o neutron-neutron (nn) pairing gap since it is not spread deeply to the interior 
region. Thereby this gap is taken the same as in our previous works, cf. |2l]. In 
[2^ |29] we have also demonstrated that the decline of the cooling curve describing 
the evolution of the Cas A surface temperature is sensitive to the value of the heat 
conductivity. In [22] we demonstrated that at an appropriate choice of the proton 
pairing gap (following model I) we are able to £t the 4% decline ACIS-S data on Cas 
A using the same lepton heat conductivity as in [30]. With the gaps given by model 
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II using the same lepton heat conductivity we match 2% decline. 

Here, we adopt the cooling inputs such as the neutrino emissivities, specific heat, 
crust properties, etc., from our earlier works performed on the basis of the HHJ EoS 
[25] and the HDD EoS |2S] for hadronic matter. The heat conductivity is the same 
as in [29]. The best £t of Cas A data with the HDD EoS was obtained in [29] with 
the same effective pion gap and the same IS'o pp pairing gap of the model I, as in our 
previous works [2H [25l |26]. Now exploiting the DD2 EoS to get the best £t of the 
cooling data we will additionally tune the pp pairing gap and the pion effective gap 
retaining all other values the same as in [29] . 



Figure 2: Square of the effective pion gap as a function of the density without pion 
condensation (curves la+lb), ttItt is the pion mass. The dotted line corresponds to 
the same parameterization as in our previous works, the solid line demonstrates a 
stronger softening effect. 

The density dependence of the square of the effective pion gap uj*'^{n) that we 
exploit in the given work is shown in Fig. [^ To be specihc we consider the case when 
the pion softening is saturated and pion condensation does not occur. The dotted 
curve la+lb is precisely the same as in Fig. 1 of [21], demonstrating saturation of 
the pion softening for n > Stiq. The solid line shows a stronger pion softening effect 
with a saturation for n > 2.5no. 

The IS'o PP pairing gaps exploited by different authors are shown in Fig. j^for the 
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Figure 3: IS'o pp gaps as functions of the proton density. The abbreviations in 
the legend correspond to those used in Ref. |36]- The gaps labeled as “Yak” and 
“AV18” are those (models I and II, respectively) exploited in our previous works 

[21I2S1I2SII2S]. 


DD2 EoS. We use the parametrization of the pp pairing gaps, Ap(pF,p), from [36], Eq. 
(2), pF,i denotes the Fermi momentum of the species i. The parameters are taken 
to fit the gaps computed in various publications. The abbreviations of the curves in 
Fig. I are taken over from Table II of [36]. Two additional curves labeled as “Yak” 
and “AV18” correspond to the models I and II, respectively, exploited in our previous 
works [211122112S] for the HHJ EoS and in [21] for the HDD EoS. 

With these gaps and the u* parametrizations we compute the NS cooling history. 


4 Results 

In Fig. we demonstrate the NS cooling history computed with the HDD EoS in 
Ref. [21] using model I for the pp pairing gaps (see the gap “Yak” in Fig. and with 
the effective pion gap given by the dotted la+lb line in Fig. As we see, the Cas 
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Iogio(t[yr]) 


Figure 4: Cooling curves for a NS sequence according to the hadronic HDD EoS; Tg 
is the redshifted surface temperature, t is the NS age. The effective pion gap is given 
by the dotted curve la+lb in Fig. The ISq pp pairing gap corresponds to model 
I. The mass range is shown in the legend. Cooling data for Cas A are explained with 
a NS of mass M = 1.497 Mq. A comparison with Cas A ACIS-S and HRC-S data is 
shown in the inset. 

A data are described by a NS with the mass McasA — 1.497Mq. A slight change 
of the value McasA compared to the value 1.54Mq found in [22] is due to inessential 
modihcations of the parametrization in the present work. 

The stiffer DD2 hadronic EoS compared to those for the HHJ and HDD EoS 
produces a smaller central density for the star of the given mass and therefore it 
leads to a weaker cooling activity, provided the same inputs are used for the effective 
pion gap u:*{n), the pairing gaps and the other model ingredients. As the result, a 
larger NS mass is required to describe Cas A cooling with the stiffer EoS. Ref. [55] 
demonstrated that when changing the EoS a description of all cooling data is possible 
even without changing any of the formerly adjusted cooling inputs except a tuning of 
the heat conductivity (in line with the strategy applied before in Ref. 1251), but the 
NS mass that is required to £t the Cas A cooling data with very stiff DD2 EoS then 
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amounts to M = 2.426 Mq. 



Iogio(t[yr]) 


Figure 5: Same as Fig. |^but for the hadronic DD2 EoS and the 15o pp pairing gap 
corresponding to the model EEHOr. The cooling data for Cas A are now explained 
with a NS of mass M = 1.950 Mq. 

Now exploiting DD2 EoS we take the heat conductivity the same as in |29] (with¬ 
out any additional tuning) but tune the effective pion gap and the pp pairing gap. 
This allows to describe Cas A cooling data by a NS with a lower mass. We performed 
calculations with all the gap curves shown in Fig. and with u*{n) given by the 
dotted and solid curves la-flb (without pion condensation) in Fig. 

The resulting cooling curves are shown in Fig. for the pp pairing gap of the 
model EEHOr and for the effective pion gap given by the dotted curve la+lb. With 
other pp pairing gaps we get a less regular description of the cooling data. The cooling 
of the hot source XMMU-J1732 is explained by a NS with the mass ~ 1.5 Mq. At 
rather low densities relevant for stars with a mass M ~ Mq the pion softening effect 
is not pronounced. A large pp proton gap at such densities is required, otherwise 
the cooling curves go down. The description of the coldest objects requires a pp gap 
dropping to zero at central densities reached in those massive objects. The cooling is 
determined by the efficient medium modihed Urea process (in the absence of pairing). 
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The description of the steep decline of the cooling curve for Cas A as indicated by the 
ACIS-S requires a small heat conductivity (in our case the lepton heat conductivity 
is computed following [SU] and the nucleon contribution incorporates a decrease with 
increasing density owing to the pion softening effect) and an efficient medium modihed 
Urea process. Cooling data for Cas A are explained with a NS of M = 1.950 Mq. 



Iogio(t[yr]) 


Figure 6: Same as Fig. but with the effective pion gap given by the solid curve 
la+lb in Fig. The cooling data for Cas A are now explained with a NS of mass 
M = 1.797 Mq. 

Fig. [^demonstrates the same as Fig. j^but with the effective pion gap calculated 
with the help of the solid curve la+lb in Fig. We see that the Cas A cooling data 
can now be explained with a lower NS mass, of M = 1.797 Mq. The cooling of the 
hot source XMMU-J1732 is also explained by a NS with a smaller mass, ~ 1.4 Mq. 
These values could be still decreased if we assumed a steeper u* {n) dependence than 
that we have chosen. Note that variational calculations of Ref. |H] show that pion 
condensation in neutron star matter may appear already for n ~ 1.3no, been in favor 
of a steeper uj*{n) dependence. Just in order to be conservative as much as possible 
we continue to exploit a much weaker pion softening in our calculations. 

In Fig. we show the same as in Fig. but for the IS'o pp pairing gaps corre- 
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Iogio(t[yr]) 

Figure 7: Same as Fig. but with the IS'o pp pairing gap corresponding to the 
model CCYms. The cooling data for Cas A are now explained with a NS of mass 
M = 1.794 Mq. 


sponding to model CCYms. With these gaps we reproduce HRC-S data for the NS 
mass M = 1.794 Mq, which is only slightly different from M = 1.797 M©, with which 
we reproduced ACIS-S data. 


5 Concluding remarks 

As we demonstrated with the DD2 EoS, large values of the NS radii and the maximum 
mass, as they might be motivated by observations, are compatible with our nuclear 
medium cooling scenario provided one uses a stiff EoS. The presently known cooling 
data on Cassiopeia A and the hot source XMMU-J1732, as well as other cooling 
data, are appropriately described by purely hadronic stars within the nuclear medium 
cooling scenario, under the assumption that different sources have different masses. 
The resulting cooling curves are sensitive to the value and the density dependence of 
the pp pairing gap and the effective pion gap. Choosing the pp pairing gap such that it 
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is large for densities relevant to low mass objects and disappearing for higher densities 
met in centres of more massive stars we are able to reach an overall agreement with 
the cooling data including Cassiopeia A and the hot source XMMU-J1732, exploiting 
soft as well as stiff hadronic EoS. Exploiting a stiff EoS and allowing the effective 
pion gap to decrease with increasing density we are able to diminish the value of the 
mass of the neutron star in Cassiopeia A required for an optimal description of its 
cooling data. Fitting a steep decline of the cooling curve for Cassiopeia A requires a 
rather low value of the heat conductivity, an appropriate form of the pp gap and an 
efficient medium modihed Urea neutrino emissivity. 
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